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Abstract 

A complete classification of exact solutions of ghost-free, massive 
bigravity is derived which enables the dynamical decoupling of the 
background, and the foreground metrics. The general decoupling so¬ 
lution space of the two metrics is constructed. Within this branch of 
the solution space the foreground metric theory becomes general rel¬ 
ativity (GR) with an additional effective cosmological constant, and 
the background metric dynamics is governed by plain GR. 

Keywords: Non-linear theories of gravity, massive gravity, bigravity, 
dark energy. 
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1 Introduction 

Recently, a ghost-free HIES nonlinear massive gravity theory was constructed 
in [31 E]. This theory is a nonlinear generalization of |5]. Independently 
later on, this ghost-free massive gravity with a flat reference metric was also 
extended to include a general background metric in [6l [7118] . A ghost-free 
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two-dynamical-metric theory, namely the bigravity as a covering theory of 
the massive gravity has also been proposed by introdncing the dynamics for 
the backgronnd metric [a [ini cn eg. 

In this paper, by referring to the simple observation already pointed 
ont in [To] which leads to a dynamical deconpling of the backgronnd [/], 
and the foregronnd [g\ metrics we will derive the general solntions f{x^) = 
F{g{x^),x^) which enable the two metrics to be solntions of two disjoint 
general relativity (GR) theories. This is possible if a portion of the effective 
energy-momentnm tensor entering into the ^f-metric eqnations as a course of 
the interaction Lagrangian between the two metrics vanishes. As the same 
term also appears in the /-metric equations by being the only contribution, 
when it vanishes the sets of held equations of the two metrics completely 
decouple from each other yielding only an algebraic matrix equation which 
generates this picture. This matrix equation written for /, and g plays the 
role of a solution ansatz that leads us to a branch of the solution space gen¬ 
erated by a Cartesian product of two GR’s. This matrix ansatz equation 
will be at the center of our analysis. In the following, we will derive the 
general solutions of this cubic matrix equation when none of the parameters 
of the theory vanish. Thus, we will be able to give a complete description of 
the solution space r[/, ^f] whose elements lead to this dynamical decoupling 
of the two metric sectors. We will also show that, the classiheation scheme 
of the analytically available solutions {f,g} admits a similarity equivalence 
class structure. 

In Section one, following a summary of the bigravity dynamics we will ob¬ 
tain the decoupling ansatz matrix equation we have mentioned above. Then, 
in the next section, we will derive the general solutions of this cubic matrix 
equation for generic constant coefficients. Since the coefficients in the actual 
equation are functions of the elementary symmetric polynomials of the so¬ 
lutions themselves a more rehning analysis is needed. Therefore, later on, 
we will present a parametric derivation which enables us to construct not 
only the solutions of this involved matrix equation, but also their elemen¬ 
tary symmetric polynomials as functions of the parameters of the theory. 
Subsequently, we will show that a subset of the generic solutions for con¬ 
stant coefficients must be omitted, when one insists on having an entire set 
of nonzero theory parameters. Besides, some of the generic solutions are 
forced to yield the same form when they are plugged into the actual matrix 
equation we have. In Section four, we will also present a formal dehnition of 
the decoupling solution space T[f,g] of the bigravity theory. We will show 
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that, this space contains a major subset that is composed of analytically well- 
dehned similarity equivalence classes of solutions. Finally, in Section hve we 
will explicitly construct the proportional background solutions, and give an 
example of Friedmann-Lemaitre-Robertson-Walker (FLRW) on FLRW case. 


2 The dynamical decoupling 


The action for the ghost-free bimetric gravity [a EQi m [12] for the fore¬ 
ground, g, and the background, / metrics in the presence of two types of 
matter can be given as 




IQttG 

K 

IQnG 


j 


R3 + A9 -2m^Cint{V^) 


+ Si 


M 


dx^ 




+ 


(O 


where R^, R-^, , A-^ are the corresponding Ricci scalars, and the cosmologi¬ 
cal constants for the two metrics, respectively. S{^ are the two different 
types of matter which independently couple to g, and /, respectively. The 
interaction Lagrangian of the two metrics is 

= Aei(\/S) -F /?2e2(\/F) + /?3e3(\/S), (2.2) 


where {e„} are the elementary symmetric polynomials 

Cl = ei(\/S) = fr\/S, 

62 = e2(\/F) = ^((fr\/S)^ - fr(\/S)^), 

63 = 63 (\/F) = \ - 3 trVf. + 2 , (ESj) 

0 

of the square-root-matrix 

= V^- (2.4) 

Likewise in [12] the original interaction terms = /Sq, and /5464 = /l 4 det\/S 
are trivially plugged into the cosmological constants A^, and Aj, respectively. 
If we demand that Eq. fl2.2l) gives the Fierz-Pauli form in the weak-held limit 
then we must have [9] 

jdi + 2 /I 2 + 1^2, = —1. (2.5) 
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Varying Eq. fl2.ip with respect to g gives the ^f-equation 

RU - Ir’S,. - = SnGTI,^. ( 2 . 6 ) 

Also, variation with respect to / resnlts in the /-equation 

K - m%i = eSnGTl,^^. (2.7) 

In these held equations the contributions coming from the interaction term 
that is given in Eq. fl2.2p are the effective energy-momentum tensors 

TJj. = gt^pK - ^intgpu, (2.8) 

and 

(2.9) 

respectively. Here are the elements of the matrix r |T2] 

r = /33(\/E)^ — (/32 + /33ei)(\/E)^ -I- (/3i + /? 2 ei -I- /33e2)\/E, (2.10) 

namely = [tY^- Both of the effective energy-momentum tensors must be 
covariantly constant 

VJ(rT„ = 0, Vl(T>r, = 0. (2.11) 

If one of these constraints is satished then the other one is automatically 
satished dSUH]. As discussed in [10] if one chooses 

r = 0, (2.12) 

then dynamically the Eqs. (12.61) . and (12.7p decouple from each other. In this 
case the hrst of the constraints fl2.1ip gives 

d^Cint = 0, (2.13) 

and thus, as the solution we will take 

C^nt = -^A. (2.14) 
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(2.15) 


Therefore, from Eqs. fl2.8p . fl2.9p we have 

' 7 ;{, = 0 . 

Consequently, the (^-equation Eq. fl2.6p becomes the usual Einstein equa¬ 
tions for g with an additional effective cosmological constant A, whereas the 
dynamically-disjoint /-equation Eq. fl2.7p reduces to be coupling-constant- 
modified Einstein equations for /. The rest of our analysis will be devoted 
to hnd the general solutions of the matrix equatioiJi] 

+ S(\/E) + ^ 14 ^ = 0, (2.16) 

where 


A = ~/^3, 

B = 132 + /^sCi, 

C = —(3i — /?2ei — ( 3362 , (|2.17|) 

which constitute the effective solution space r[ 5 f, /] of the ghost-free bigravity 
action (12. ip that enables the above-mentioned dynamical decoupling for the 
foreground, and the background metrics. 

3 The structure of the solution space P 

Now, let us consider the matrix equation 

+ BX'^ + CX = X{AX^ + BX + ^ 14 ) = X{X - Ail 4 )(X - A 2 I 4 ) = 0, 

(3.1) 

for a 4 X 4 matrix function X{x^). For the following analysis we will disre¬ 
gard the solutions which require either of the /3-coefficients to be zero. The 
characteristic polynomial of any 4x4 matrix X would be the degree-four 
polynomial 

Pxit) = detitli — X), (3.2) 

whose four roots are the eigenvalues of X. We should observe first that if X is 
a solution of Eq. fl3.ip then for any invertible matrix function P{x^), P~^XP 

^We will prefer to work with the negative of the Ea. (l2.12p in accordance with the 
massive gravity formalism [mile], and for future relevance. 


5 















is also a solution. Therefore, in order to find the general solutions of Eq. fl3.ip 
it would be sufficient to classify the Jordan canonical forms satisfying flJ.ip . 
Our main objective next, will be the classihcation of the similarity equiva¬ 
lence classes of solutions with respect to their minimum polynomials m{X). 
The roots of m{X) are the same with the eigenvalues of the various equiv¬ 
alence classes of matrices having that minimum polynomial with differing 
multiplicities of course. Since Eq. fl3.ip is a degree-three polynomial equation 
when it is factorized, its various factors with degrees smaller than or equal 
to three will dehne the minimum polynomials of its 4 x 4 matrix function 
solutions. In other words, the solutions can be classihed with respect to these 
minimum polynomial factors. 

3.1 The algebraic structure 

In the following classihcation, we will identify the entire set of similarity 
equivalence classes of solutions satisfying Eq. fl3.ip by simply taking the coef- 
hcients in Eq. fl3.1l) to be constants. We will group the solutions with respect 
to their minimum polynomials and it will be sufficient to determine the Jor¬ 
dan canonical form spectrum of each minimum polynomial which corresponds 
to some combination of the factors in Eq. fl3.ip . 

1. Solutions with m{X) = X 

Since now, the minimum polynomial has no repeated roots and its 
unique root is zero these solutions are diagonalizable with zero eigen¬ 
values. They are the trivial solutions X = 0. These solutions would 
demand at least one of the metrics to be zero via Eq fl2.4l) as metrics 
are invertible matrices. 

Next, we will classify the solutions of Eq. fl3.ip with respect to the root 
structure of the quadratic factor. 

• The cases when — AAC > 0 : 

In this case, the factor AX‘^ -|- BX + CI 4 has two distinct real roots 
-^ 2 - 

2. Solutions with m{X) = X — Ai^ 2 l 4 

Now, there exist diagonal solutions as the minimum polynomials have 
no repeated roots and they are linear. These solutions are 

Ui = A 1 I 4 , U 2 = A 2 I 4 , (3-3) 
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with minimum polynomials m{X) = X — A1I4, m{X) = X — A2I4, 
respectively. 


3 . Solutions with m{X) = AX"^ + BX + 

The matrices with a minimum polynomial of the form m{X) = AX'^ + 
BX + 6*14 are diagonalizable when the roots are distinct like our case 
here. Thus, by considering all the possible multiplicities of the eigen¬ 
values which are the same with these distinct roots Ai,A 2 the Jordan 
forms corresponding to this minimum polynomial become 



/Ai 

0 

0 

0\ 


/Ai 

0 

0 

0\ 

U 3 = 

0 

A 2 

0 

0 

,U,= 

0 

Ai 

0 

0 

0 

0 

Ai 

0 

0 

0 

Ai 

0 


VO 

0 

0 

1 ^ 2 / 


Vo 

0 

0 

A 2 / 


/A 2 0 0 0\ 

0 A 2 0 0 

0 0 A 2 0 

\0 0 0 Ai/ 


dsaD 


4 . Solutions with m(X) = X{X — A1I4) 

In this case, the roots of the minimum polynomial are {0, Ai}. Again, 
the matrices with this minimum polynomial are diagonalizable. The 
possible Jordan forms with {0, Ai} eigenvalues read 


Ue = 


/Ai 0 0 0\ 

0 0 0 0 
0 0 0 0 
\0 0 0 0 / 




Ai 
0 
0 
VO 


0 

Ai 

0 

0 


0 o\ 
0 0 
0 0 
0 0 / 


,Us 


/Ai 

0 

0 

0\ 

0 

Ai 

0 

0 

0 

0 

Ai 

0 

^0 

0 

0 

0/ 




(3.5) 


5. Solutions with m(X) = X{X — A 2 I 4 ) 

Similarly, the diagonalizable solutions with eigenvalues {0, A 2 } as the 
roots of the minimum polynomial have the following possible Jordan 
forms 



( A 2 

0 

0 

0\ 


( A 2 

0 

0 

0\ 


( A 2 

0 

0 

0\ 

Ug = 

0 

0 

0 

0 

,Uio = 

0 

A 2 

0 

0 


0 

A 2 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

A 2 

0 


Vo 

0 

0 



Vo 

0 

0 

V 


Vo 

0 

0 

0/ 














(3 

.6) 
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6 . Solutions with m{X) = X{X — Ail 4 )(X — A 2 I 4 ) 

Now, the roots of the minimum polynomial are {0, Ai,A 2 }. They are 
distinct, and this states that the solutions whose minimum polynomial 
is of this form must be again diagonalizable. The possible Jordan forms 
constructed from these eigenvalues are 


Ul2 = 


/o 

0 

0 

Vo 


0 0 
Ai 0 
0 Ai 
0 0 


17i4 = 


0\ 


/o 

0 

0 

0\ 

0 


0 

0 

0 

0 

0 

0 

0 

Ai 

0 

X2J 


VO 

0 

0 

X2/ 

/o 0 0 0 \ 





0 Aa 0 0 

0 0 Aa 0 

\0 0 0 Ai/ 




• The cases when — AAC = 0 : 




For these cases, the factor AX"^ + BX + CI 4 has a repeated real root 
which we will call A'. 


7. Solutions with m{X) = X — A '14 

This solution, is another diagonal one due to its minimum polynomial. 
It is 

Vo = X'u. (3.8) 

8 . Solutions with m(X) = (X — A' 14 )^ 

Since now, the minimum polynomial has repeated roots the matrices 
with such a minimum polynomial are nondiagonalizable. Their eigen¬ 
values must be A' with multiplicity four. Thus the primary block is four 
dimensional and since the multiplicity of the roots of m{X) is two the 
maximum dimension of the secondary block must be two. The possible 
Jordan forms are 











/A' 

1 

0 

o\ 


/A' 

0 

0 

o\ 

0 

A' 

0 

0 

, J"2 = 

0 

A' 

1 

0 

0 

0 

A' 

0 

0 

0 

A' 

0 

VO 

0 

0 

^7 


VO 

0 

0 

^7 

/A' 

0 

0 

0\ 


/A' 

1 

0 

0\ 

0 

A' 

0 

0 

, 14 = 

0 

A' 

0 

0 

0 

0 

A' 

1 

0 

0 

A' 

1 

^0 

0 

0 

^7 


Vo 

0 

0 

^7 



dSSl) 


We should remark that, it would be enough to take either of {Vi, V 2 , V 3 } 
however, to be on the safe side we include all of them. 

9. Solutions with m{X) = X{X — A' 14 )^ 

The roots of the minimum polynomial in this class of solutions are 
{0, A', A'}. Once more, m{X) has repeated roots so the matrices with 
this minimum polynomial are nondiagonalizable. The dimension of 
the secondary blocks for the eigenvalues A', and 0 are two, and one, 
respectively. The possible Jordan canonical forms are 


/A' 

1 

0 

0\ 



/A' 

0 

0 

0\ 

0 

A' 

0 

0 

,^6 


0 

A' 

1 

0 

0 

0 

A' 

0 

— 

0 

0 

A' 

0 

Vo 

0 

0 

V 



Vo 

0 

0 

0/ 




/A' 1 

0 

0\ 





^7 


0 A' 

0 

0 





— 

0 0 

0 

0 







\0 0 

0 

V 





Like the previous case, we include both V 5 , and Vq in spite of the fact 
that one would be enough to generate the similarity class of solutions. 

10. Solutions with m{X) = X{X — A' 14 ) 

There are no repeated roots of m{X) for these solutions. Therefore, 
the matrices with this minimum polynomial are diagonalizable. The 
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possible Jordan normal forms with the eigenvalues {0, A'} are 


/A' 

0 

0 

0\ 



/A' 

0 

0 

0\ 

0 

0 

0 

0 

,^9 


0 

A' 

0 

0 

0 

0 

0 

0 

— 

0 

0 

0 

0 

Vo 

0 

0 

V 



Vo 

0 

0 

0/ 




/A' 

0 

0 

0\ 







0 

A' 

0 

0 






0 

0 

A' 

0 







VO 

0 

0 

V 





• The cases when — AAC < 0 : 




For this class of solutions, the factor AX"^ + BX + has complex 
roots which we will call A, A*. In this respect for such cases, the mini¬ 
mum polynomial of the similarity equivalence classes of real solutions 
must include the quadratic factor as a whole since for real matrices 
if A is an eigenvalue A* must also be an eigenvalue, and vice versa. 
For this reason, there are two possible minimum polynomials and the 
corresponding solutions are nondiagonalizable. 


11. Solutions with m{X) = {X — Al 4 )(X — A*l 4 ) 


The eigenvalues of the matrices with this m{X) are {A, A*}. If we dehne 
the real and the imaginary parts of A = then the unique Jordan 

canonical form for the real solutions with this minimum polynomial 
becomes 


/ R I 0 0\ 

-I R 0 0 

^ ~ 0 0 R I ' 

\o 0 -I rJ 

12. Solutions with m(X) = X{X — Al 4 )(X — A*l 4 ) 


(3.12) 


The eigenvalues of the solutions of this equivalence class are {0, A, A*}. 
The only possible Jordan normal form for the real matrices satisfying 
m{X) = X{X — Al 4 )(X — A*l 4 ) = 0 is of the form 


/o 0 0 o\ 
0 0 0 0 
0 0 R I 
\0 0 -I Rf 


(3.13) 
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For the purpose of achieving the most general construction of the so¬ 
lution space we should also consider the cases when the coefficients 
in Eq. fl3.ip vanish. We will disregard the cases when A = 0 as they 
correspond to choosing /Sa = 0 . 

• The cases when i? = 0 : 


If we assume that CjA < o[§ then the quadratic factor in Eq. fl3.ip can 
be written as {X — A)(X -|- A) where A = ^J—CjA. 

13. Solutions with m{X) = X ± AI 4 

There exist diagonal solutions. They are 


Ki — AI 4 , K 2 — —AI 4 , (3.14) 


with the minimum polynomials m{X) = X— AI 4 , and m(X) = X-I-AI 4 , 
respectively. 


14. Solutions with m{X) = AX^ + CI4 

The roots of the minimum polynomial thus, the eigenvalues in this case 
are {A, —A}. They are distinct hence, the matrices with this minimum 
polynomial are diagonalizable. The possible Jordan forms are 


/A 

0 

0 

0 \ 


A 

0 

0 

0 \ 

0 

A 

0 

0 

,/^4 = 

0 

A 

0 

0 

0 

0 

A 

0 

0 

0 

-A 

0 

\0 

0 

0 

-A/ 


\0 

0 

0 

-A/ 


/AO 0 0 \ 

0 -A 0 _ 0 

0 0 -A 0 

VO 0 0 -A/ 


15. Solutions with m(X) = X(X — AI 4 ) 


(I3T5]) 


The roots of the minimum polynomial and the eigenvalues of the so¬ 
lutions in this case are {0,A}. They are not repeated so, the cor¬ 
responding matrices are diagonalizable. The Jordan canonical forms 


^We will comment on the case when C/A > 0 in the next subsection. 
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with various eigenvalue multiplicities are 


/A 

0 

0 

o\ 


/A 

0 

0 

0\ 

0 

A 

0 

0 

,^7 = 

0 

A 

0 

0 

0 

0 

A 

0 

0 

0 

0 

0 

\0 

0 

0 

0/ 


\o 

0 

0 

0/ 


/A 0 0 0 \ 
0 0 0 0 
0 0 0 0 
\0 0 0 0 / 


(EIED 


16. Solutions with m{X) = X{X + AI 4 ) 


Now, the roots of the minimum polynomial and the eigenvalues of the 
solutions become {0, —A}. Again, they are distinct, the corresponding 
matrices are diagonalizable. The possible Jordan canonical forms with 
various eigenvalue multiplicities are 


/-A 

0 

0 

0\ 


/-A 

0 

0 

o\ 

0 

-A 

0 

0 


0 

-A 

0 

0 

0 

0 

—A 

0 

0 

0 

0 

0 

V 0 

0 

0 

0/ 


\ 0 

0 

0 

0/ 


/-A 0 0 0 \ 
0 0 0 0 


\ 0 0 0 0 / 


(EtZI) 


17. Solutions with m(X) = X{AX'^ + CI 4 ) 

In this case, the roots of the minimum polynomial are {0, A, —A}. They 
are not repeated, hence, the corresponding matrices are diagonalizable. 
The possible Jordan canonical forms constructed with various multi- 
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plicities of these eigenvalues can be listed as 


Ki2 = 


(X 

0 

0 

o\ 


(X 

0 

0 

0 \ 

0 

-A 

0 

0 


0 

A 

0 

0 

0 

0 

0 

0 

0 

0 

—A 

0 

Vo 

0 

0 

0 / 


Vo 

0 

0 

0 / 


Ku = 


(X 0_ 
0 -A 
0 0 
Vo 0 


0 o\ 
0_ 0 
-A 0 
0 0 / 


(M) 


The cases when C = 0 : 


18. Solutions with m{X) = AX + 5 I 4 

Since, the minimum polynomial is linear this case is a diagonal solution. 
It is 


Xi — (3.19) 

19. Solutions with m{X) = X[AX + 5 I 4 ) 

The roots of m{X) now, become {0, —5/^4}. They are distinct thus, 
the matrices with this minimum polynomial are diagonalizable with 
eigenvalues {0, —B/A}. The possible Jordan normal forms are 


X2 


(-B/A 

0 

0 

0\ 


(-B/A 

0 

0 

o\ 

0 

0 

0 

0 

0 

0 

0 

0 

,^3 = 

0 

0 


-B/A 

0 

0 

0 

0 

0 

V 0 

0 

0 

V 


V 0 


0 

0 

0/ 


/- 

B/A 

0 

0 

o\ 




Y — 


0 


-B/A 

0 

0 




VV4 — 


0 


0 

-B/A 

0 





V 

0 


0 

0 

V 





(E2ni) 


20 . Solutions with m(X) = X‘^[AX + BI 4 ) 
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For this branch of solutions, the minimum polynomial m{X) has re¬ 
peated roots which are {0, 0, —S/A}. Hence, the matrices with this 
minimum polynomial are nondiagonalizable. In the relevant Jordan 
canonical forms, the dimension of the secondary block corresponding 
to the zero eigenvalue is two, and for the —B/A eigenvalue it is one. 
The possible Jordan forms are 


^5 


/-B/A 

0 

0 

o\ 



/- 

-B/A 

0 

0 

0 \ 

0 

0 

1 

0 

,Xe = 



0 

0 

0 

0 

0 

0 

0 

0 



0 

0 

0 

1 

V 0 

0 

0 

oj 



V 

0 

0 

0 

0 / 


/- 

B/A 

0 

0 

o\ 





^7 = 


0 

0 


-B/A 

0 

0 

0 

0 

1 






V 

0 


0 

0 

V 
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3.2 The parametric structure 


mu 


In the previous subsection, we have derived and classihed all the repre¬ 
sentatives of the similarity equivalence classes of the solutions of Eq. fl3.ip . 
The reader may verify that the matrices listed for various conditions satisfy 
Eq. fl3.ip by direct substitution for unspecihed constant coefficients under 
these conditions. Therefore, all the matrix functions constructed from these 
representatives by similarity transformations via an invertible matrix func¬ 
tion P{x^) also satisfy Eq. fl3.ip . However, in Eq. fl2.16p the coefficients of the 
matrix equation for \/E are not numerical constants, on the contrary, they 
are functions of the elementary symmetric polynomials 61,62 of \/E which 
turn this matrix equation into a very nontrivial one. In this respect, our 
task is not over and we have to furthermore solve the entries of the matrices 
listed in the previous subsection explicitly, so that they satisfy this nontrivial 
matrix equation Eq. fl2.16p . 


• The solutions; t/1,2,3,4,5, To,i,2,3,4> ^ = 


In this case, the solutions satisfy 

AX^ + BX + CI4 = 0. 


(3.22) 
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By taking the trace of this equation, and using Eqs. fl2.3p . and fl2.17p we get 


3/^261 + 2/^362 + 4/?i — 0. (3.23) 

For the cases which satisfy — 4AC > 0 the roots of Eq. fl3.22p are 
B + y/W^^4AC 


Ai = 

From Eqs 


Ao — 


-B - VB^ - 4AC 


2 A ’ 2 A 

, and fl3.4p we have 

Cl = tr{Ui) = uAi + mA 2 , 


(3.24) 


(3.25) 


for f/i, t/ 2 , f/4, f/5 the (u, m) values are (4, 0), (0,4), (2, 2), (3,1), (1, 3), re¬ 

spectively. When m — n 7 ^ ±2 namely, for Ui,U 2 , U 3 using Eqs. (13. 23 p . and 
(I3.24P in Eq. (l3.25p yields the quadratic equation 


o(ei) -|- bci "he — 0, 


(3.26) 


with 


a = 2(n - l)(m - l)(/93)^ 
b = — ((m — n)^ -f 2 (m + n — 2 mn))(] 2 l 33 , 


c = 2 mn{f 32 Y — 2 {n — mYf 3 if 33 . 


Therefore, provided 


— 4ac > 0, 

(3.28) 

Cl becomes 


—b ± — 4ac 

(3.29) 

“ 2 a 


and from Eq. (13.231) we have 


3/^2 -4ac^ 2/3i 

- 2 ;^-) - w 

For U 3 , 5^ — 4ac = 0 thus, the condition fl3.28p is automatically satished. For 
Ui, and U 2 it becomes 

3{(32Y - 4 / 3 i /33 > 0 . ( 3 . 31 ) 
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Furthermore, we can now explicitly write the components of Ui,U 2 , as 


+ . _/32 + /33(-^±^f^)+\/A 

% ’ 

(3.32) 

where 

A ± \/62 — 4acx2 , ^ ^ /—b± y/b'^ — 4acx _ ,n ^ ^ ^ 

“ “ (3.33) 

In this expression, by substituting a, b, c, from Eq fl3.27p the restriction A = 
— AAC > 0 that dehnes these solutions takes the form 

-3(/32)' + 4A/93>0, (3.34) 

for f/ 3 , and it reads 

3(/32)'-4/3i/33>0, (3.35) 

for Ui, and 1 / 2 - The condition fl3.35p is stronger than fl3.3ip hence, it must 
be taken as the dehning constraint on the parameter space for the solutions 
Ui,U 2 to exist. Therefore, when the coefficients {/?*} satisfy Eg. (13.341) the 
similarity equivalence class of (upon the substitution of the matrix ele¬ 
ments via fl3.32lF are solutions of Eq. fl2.16p . Also, when the coefficients {/S*} 
satisfy Eq. fl3.35p the similarity equivalence classes of Ui, and U 2 again, by 
using the relevant matrix entries from fl3.32p are solutions of Eq. fl2.16p . On 
the other hand, when m — n = ±2, namely, for the cases t /5 Eqs. fl3.23p . 
fl3.24p . and fl3.25p lead us to the condition 

3(/?2)'-4/3i/33 = 0. (3.36) 

By using this condition one can show that we always have 

A „±2 = 52 - 4AC = (2/^2 + > 0. (3.37) 

Thus, for these cases A > 0 is always satisfied provided fl3.36p holds. The 
roots of the minimum polynomial now become 

, _ /^2 A (A^ei — (±( 2/?2 + /^s^i)) . _ 1^2 + /^sCi -f (±(2/32 ± PzA)) 

- % ’ - % ■ 

(3.38) 
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By using these in t/ 4 , and t /5 the consistency condition 03.251) demands that 
the signs in 03.38p must be chosen such that 


f /4 


U, 


/ P 2 
2/33 


0 

^ 2/33 


0 

0 


/M 


0 
0 

'2+2/3361 


0 

0 


02 


0 

0 

0 




203 

0 

0 

0 


2/^3 , , 

3fe+2^3ei 
^ 203 / 

0 0 0 \ 
0 


0 

^ 2/33 


0 


0 


2 / 33 / 


(iraD 


which are similar to each other thus, we will take the representative of the 
equivalence class as U = U 4 . We remark that, in this case the formulation 
does not specify Ci and it remains as a free parameter. It can take any value 
except — 2 ( 32/(33 which would violate 03.37p . Therefore, the solution U = f/ 4 , 
is parametrized by a free trace parameter ei as in Eq. 03.39p . and from 03.23P 
we have 


3/^2 2A 
/Ss ■ 


(3.40) 


Among the solutions which satisfy Eq. 03.22p : when — AAC = 0 we have 
the ones Vq, fh, V 2 , V 3 , V 4 , and when B'^ — AAC < 0 we have the solution Y. 
For all of these cases taking the trace of these matrices leads us to 


ei = 


/Ss ■ 


(3.41) 


From Eq. (l3.23p we get 

^{^ 2 ? 2(3, 

{h? ft ■ 


(3.42) 


From these relations for Vq, Vi, V 2 , V 3 ,14 via the dehnitions in Eq (l2.17p we 
have 


2 A 2/33’ 


(3.43) 


and the condition B‘^—AAC = 0 which allows these solutions to exist becomes 


3(42)' - 44 i 43 = 0. 


(3.44) 
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On the other hand, when — AAC < 0 is valid Eq. fl3.22p has two complex 
roots X = R + li, and \* = R — li, here 


1^2 j _ \/3(/?2)^ — 4/3i/^3 

Wz ~ 


(3.45) 


where we have used Eqs. 02.171) . 03.4ip . and 03.42p . From the same equations 
the condition — AAC < 0 which allows the solution Y to exist now, 
becomes 

-3(/32)' + 4A/d3<0. (3.46) 

• The solutions; hs ,,, lo, • 

For the 7/—matrices 


ei = nXi + mX2 


— {n + m) 


B 


+ (n 


m) 


V52 - AAC 
2A 


(3.47) 


where (n,m) = (1, 0), (2, 0), (3, 0), (0,1), (0, 2), (0, 3), (2,1), (1,1), (1, 2) for 
f/g, f/y, f/g! Ug, f/io, f/ii, f 7 i 2 , h^i 3 , h^i 4 , respectively. We also have 

e 2 = ]^{{eif -tr{Uif). (3.48) 

Now, from this relation by using tr{UiY — ''^(•^i)^+^(-^ 2 )^, and the fact that 
Ai, A 2 are the roots of Ax^ + Bx + C = 0 we obtain 

= (2 + <" + ’">1^’ 

when n + m ^ 2. Again, for the cases n + m 7 ^ 2 substituting Eqs. 02.17p . 
and (I3.49P into Eq. (l3.47p after some algebra leads us to the equation 

a'(ei)^ + 6 'ei + c' = 0, (3.50) 


where 


a' = (2 — u — m)(l — m)(l — n){fAzf, 

h' = —((2 — n — m)(l — m)n + (2 — n — m)(l — n)m — {n — 'm)^)/52/d3, 
d = mn{2 — n — m){l32Y A 2{n — mY fl3.5ip 
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If the condition 


72 - 4a'c' > 0, 

(3.52) 

is satished then ei becomes 


—b' ± \/&'2 — 4a'c' 

(3.53) 

“ 2a' 


and via Eq. fl3.49p we get 


1 n , ,,/?2^-&'±\/&'2-4a'c'^ /3i^ 

62 = 77 ;-+ —- ) + i'ri + 'm)—), ( 3 . 54 ) 


(2 — n — m) 


733 


^3' 


which is now written solely in terms of the parameters {/9j}’s. We can also 
write the components of the [/-matrices as 


-^ 1,2 — 


2^3 


(o „ ± \/n^-4ah/x X 

-(A + A(-- ))±'^). 


(3.55) 


where 


A = - 4AC' 


— ( ^2 + 73 ( 
1 


-V ±^/h‘ 


n 


4a'ch 


2 a' 


) ) + 473 


,-h' ± \/b'‘^ - Aa’d. 

-27- > 


(2 — n — m) 


, 4= \/['^ — 4a'c'X 

(n -h m - 1)72(- — -) (n + m) 7 ] 


( 1336 ]) 


For the existence of these solntions we mnst have A > 0, together with the 
condition fl3.52p . In spite of the fact that, Eq. fl3.50p is generally valid for 
all of the matrices; f/e, 8 , 9 , 11 , 12,14 we mnst be cantious since for some of these 
solntions it rednces to a rednndant form, or it is trivially satished. We will 
consider snch cases one by one. In particnlar, for Uq, and Ug (I3.52p is satished 
directly as for these cases 6'^ — 4a'c' = 0. However, for these cases Eg. (13. 50 p 
holds only if 7i73 = 0) thus, we will exclude the solutions Uq, Ug. For Ug, Un 
fl3.52p becomes ( 72 )^ — 7i73 A 0, and A > 0 leads to 

(72)'-7i73>0, (3.57) 


which is stronger than fl3.52p . Therefore, for the existence of the solutions 
Us,Uii, Eq. fi 3 . 57 p is the only condition. For the matrices UuiUi^ we have 
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6 '^ — Aa'c' = 0 hence, Eq. fl3.52p is satisfied. Also, now Eq. fl3.5Up turns into 
the condition 

-(/52)' + /3i/53 = 0, (3.58) 

and by using this condition for these cases, we find that 

A = (3/32 + /53ei)2 >0. (3.59) 


We observe that, for t/ 12 , if fl3.58p is satisfied we always have A > 0, this 
fact makes Eq. fl3.58p the unique condition for the existence of these solutions. 
For the cases 17i2, the roots of the minimum polynomial are 


Ai = 


/^2 + /^sCl — (±(3/?2 + /^sCl)) 


A 2 — 


h + + (A(3/32 + /^sfii)) 


(3.60) 

When we substitute them in the matrices f/ 12 , and f /14 the consistency con¬ 
dition (I3.47P constrains us to choose the signs in Eq. (l3.6np in such a way 
that 


[/12 = 17i4 = U' = 


'14 


/o 

0 

0 


0 

A 

0 

0 


0 

0 

0 


0 
0 
0 

2/32+feei 

P3 


\ 


(3.61) 


where ei remains to be a free trace parameter, as the matrix fl3.6ip trivially 
satisfies Eq. fl3.47p . It can take any value except — 3 P 2 /P 3 which would violate 
fl3.59p . For the solution U', which is parametrized by ei from Eq. fl3.49p we 
have 

2/^2 3 / 3 i 

■ 77- 

On the other hand, if we turn our attention to the cases when n -|- m = 2, 
namely, the cases; Uj, f/io, t/ia again, via tr{UiY = n(Ai)^ -|- m(A 2 )^, and the 
fact that Ai, A 2 satisfy Ax^ + Bx -|- O = 0, Eq. fl3.48p gives 


ei = 


2/di 
' (^ 2 ' 


(3.63) 


For f/i 3 , we have ei = Ai -|- A 2 however, this relation together with Eq. 
result in the constraint /32 = 0. From Eq. (13.631) in this case we must also 
have /5i = 0. Therefore, we will disregard Ui^. For the matrices f/y, f/ 10 . 


Cl — 2 Ai_2 — 2 




2A 




(3.64) 
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where + corresponds to t/ 7 , and — to t/ 10 , respectively. By referring to the 
definitions in Eq. fl2.17l) . and by using Eq. fl3.63p in Eq. fl3.64p we find that 


(A)^ 


(3.65) 


For the solutions f/ 7 , f/io, upon the substitution of ei, 62 from Eqs. (l3.63p . and 
fl3.65p the dehning condition — AAC > 0 becomes 


(A)^ > 0 , 


(3.66) 


which is automatically satished. For these cases, by referring to the Eqs. 
fl2.17p . 03.631) . and 03.65P the computation of the matrix entries reads 


A 


1,2 — 


(A)^ 


2/ji/33 T (±(A)^) 
2AA 


(3.67) 


By substituting these in the matrices 6 ^ 7 , and t/ 10 , the consistency condition 
03.63P denotes that the signs in 03.671) must be chosen such that 




Uj = U,o = U” = 


0 0 0 \ 


0 0 0 

P 2 


0 
V 0 


0 0 0 
000/ 


Let us now focus on the solutions V5,6,7,8,9,10- They have 

(ei)^ 


ei = trVi = nX', tr{Vi) = n{X!) = 


n 


(3.68) 


(3.69) 


where for V 5 , Vg, Viq n = 3, for Vj,Vq n = 2, and for Vg n = 1. Furthermore, 


62 = ^((ei)^ - tr{Vi}'^) = (3-70) 


Referring to the definitions of the coefficients in Eq. 02.17p once more, explic¬ 
itly we have 


B _ A + A^i 
^ “ 2 A 


(3.71) 


Using this relation in Eq. 03.69p gives 


^A 

(2 - n)A’ 


(3.72) 
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when n 7 ^ 2. If we consider the matrices Vj, Vg for which n = 2 this compu¬ 
tation leads to /32 = 0, hence, we also disregard these solutions. Next, substi¬ 
tuting Eq. fl3.72p in Eq. fl3.70p we also obtain 62 in terms of the /dj-parameters 
as 


n(n - 1 )(/ 32 )^ 

2(2-n)2(/33)2- 


(3.73) 


Finally, by substituting Eq. (l3.72p into Eq. (l3.7ip we can also calculate explic¬ 
itly the matrix entries as 


A' = 

(2-n)/?3- 


(3.74) 


The condition — 4AC = 0, needed for the existence of the solutions 
14 , 6 , 7 , 8 , 9,10 reads 


1 

(2 - n )2 


n 

2 — n 


n{n — 1) 
2 {n - 2)2 


(132?-m 


0 , 


(3.75) 


where we have used Eqs. (12.171) . fl3.72p . and (13.731) . In particular, for the 
case Vg; n = 1, this condition gives /3i/33 = 0 thus, we also disregard Vg. For 
the remaining cases 14 ^ 6 , 10 ) since n = 3 Eq. (l3.75p takes the form 


(/92)'-/5i/d3 = 0. 


(3.76) 


Our final case in this set of solutions is Y' which exists when — 4AC < 0. 

Starting from its definition in Eq. (l3.13p we find that 


ei = trY' = 2R, (3.77) 

where 

(3.78) 

By substituting the dehnitions of A,B from Eq. (l2.17p into Eq. (l3.77p we get 
42 = 0. Therefore, Y' will also be excluded from our solution space which 
will be constructed for nonzero coefficients in Eq. (l2.ip . 

• The solutions; Xi... 7, 77i,...p4 : 

In general, we will exclude the cases when A = 0 as they correspond to /^s = 0. 
The cases when B = 0 with C/A>Q may have the minimum polynomials 
either m{X) = AX^ + CI 4 , or m(X) = X{AX‘^ + CI 4 ). However, both cases 
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are restricted to /92 = 0. Thus, we will exclude them. On the other hand, 
when S = 0, CjA < 0 we have the solutions Here 


Cl = tr{Kj) = nX = n\/—C / A, (3.79) 

and 

~ HI 

tr{K,f = m{Xf = -(el)^ (3.80) 

where (n,m) =(4,4), (-4,4), (2,4), (0,4), (-2,4), (3,3), (2,2), (1,1), (-3,3), 
(-2, 2), (-1,1), (0, 2), (1, 3), (-1, 3) for i7i,2,...,i3,i4, respectively. For K 4 , K 12 , 
Cl = 0, from the condition B = 0 we observe that for these cases to be so¬ 
lutions we must have /92 = 0 thus, we exclude these solutions too. For the 
other cases from the above relations we have 


62 = l((ei)^ - tr(Kif) = (3.81) 


Now, if we use Eqs. fl2.17p . and fl3.8ip in Eq. fl3.79p we obtain the quadratic 
equation 

+ + + = (3.82) 

When rA — m ^ —2 this equation has the solutions 


ei = 


I 

(33 



2(2 + n2 -m)^ 


2 -f 77,2 — m 


(3.83) 


We observe that if 

Cl is real and via the relation Eq. fl3.79p the condition C/H < 0 is satished. 
For these solutions to exist we must also satisfy the condition 5 = 0, which 
reads 

[rA - (m - 2)^] {[ 52 f - 2(2 + rA - m)rAl 3 i( 33 , = 0, (3.85) 

where we have made use of the Eqs. fl2.17p . and fl3.83p . For the solutions 
7^8, Kii (13.8511 gives /di/ds = 0, which leads us to exclude these solutions. A 
close inspection denotes that when (13.8511 is satished (13.8411 is automatically 
satished, thus, Eq. fl3.85p is a stronger condition to be considered uniquely. 
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Therefore, the condition fl3.85p is the unique dehning condition for the ex¬ 
istence of the solutions Ki, K 2 , K^, K^, Kq, K^, Kq, Kiq. On the other hand, 
when — m = —2, namely, for the cases K 13 , Kn Eq. fl3.82p leads us to the 
result 



(3.86) 


02 must again be read from Eq. (l3.8ip . In this case, since A must be positive, 
via Eq. fl3.79p we see that; when ^2 > 0, the solution is i^i 4 , and when 
A//S 2 < 0, the solution is 77i3. Besides, since C/A = —n“^(ei)^ = — (ei)^ the 
condition C/A < 0 is automatically satished for these solutions. The other 
restriction B = (32 + (3zei = 0 becomes 


(/92)'-/9i/93 = 0, 


(3.87) 


which is left as the unique dehning condition for 77 i 3 , 77 i 4 to exist. For all 
these solutions discussed above, the matrix entries can be found as A = 
= 01 /n by substituting the relevant 01 , and n values case by case. 
We should state however, that when one computes 01 from Eq. fl3.83p . then 
substitutes the results in Eqs. fl3.14p . fl3.16p . fl3.17p for the cases Ki, K 2 , Kq, 
Ky, Kg, and Kig, respectively one hnds that Ki = K 2 , Kq = Kg, and K^ = 
Kiq. For this reason, we will dehne the matrices K = Ki = K 2 , K' = Kq = 
Kg, and K" = K^ = Kiq for the rest of our analysis. Now, we will turn our 
attention to the condition (7 = 0 which generates the X—series of solutions. 
Firstly we should note that when (7 = 0, in the previous subsection, the 
solutions coming from the minimum polynomials m{X) = X, and m{X) = 
X^ are excluded as they lead to the restriction (3i = 0. Besides, the solutions 
with m{X) = AX^ when B = C = 0 are also disregarded as they are 
restricted to the cases (3i = (32 = 0. For Xi we have 01 = trXi = —AB/A 
which gives 


4/?2 


(3.88) 


We also get 

62 = ^((ei)^ - tr{Xif) = ^(ei)^. (3.89) 

By using these results, and by referring to Eq. fl2.17p once more, the existence 
condition (7 = 0 becomes 


2{(32f -3>(3x(3q = Q. 


(3.90) 
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(3.91) 


Then, by making use of Eq. fl3.88p we can also explicitly write Xi as 


Xi = 


A 

3/53 


I 4 . 


For the other cases; X 2 , 3 , 4 , 5 , 6 , 7 , ei = trXi = —nB/A with n = 1 for X 2 , X 5 , Xq, 
and n = 2 for X 3 ,Xj, and n = 3 for X 4 , respectively. Thus, for these cases 
when n 7 ^ 1 we obtain 


ei = 


n /?2 


(n 




(3.92) 


For n = 1, substitution of A,B via Fq. (l2.17p in ei = —nB/A gives (32 = 0, 
hence, we exclude the solutions X 2 , X^, Xq. For the solutions ^ 3^4 7 , we can 
compute 62 as 


62 = - tr{Xif) = 

Similarly, for these cases the condition C = 0 reads 

n{/32f - 2{n - = 0. 


(3.93) 

(3.94) 


Also, by using Fqs. fl2.17p . and fl3.92p we can explicitly compute the matrix 
entries of X 3 ^ 4^7 as 


B P 2 

A (n-l)/? 3 ' 


(3.95) 


4 The decoupling solution space 


In the previous section, we have explicitly worked out and constructed the 
Jordan canonical forms of the nontrivial fixed-eigenvalue solutions of Eq. (13.ip 
in terms of the /5j—coefficients when neither of /dig ,3 is chosen to be vanishing. 
Let us first define the set composed of these solutions 


J={[A, U 2 , 1 / 3 , Us, Un, U, U', U”, Fo, F, F, F, F, F, F, Fo, F ^1, 


X 3 , X 4 , X 7 , K, K', K", K 3 , X 5 , Xi 3 , Xi 4 


(031) 


In the Appendix, we give a summary of the exact forms of these matrices 
whose entries are derived in terms of the /d—coefficients in Section three. 
Next, we define the set of matrices 

M = \mj = M-^JM \yj ej, and detM 7 ^ o|, (4.2) 
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where M is any invertible real constant matrix. As we discnssed before, 
the similarity eqnivalence classes of the elements of are also solntions of 
Eq. fl2.16p . Bearing this in mind now, the set of matrix fnnctions 

5 = |\/S -U^M I \/S e C'°°|, (4.3) 

where W is a coordinate chart on the AD spacetime manifold, is the complete 
set of general nontrivial solntions of Eq. fl2.16p when /dig,3 7^ 0. Of course, 
when the /dj-parameters are determined the appropriate elements of J must 
be chosen in other words, the ones whose domain of validity is satished by 
these parameters. We can also define the subset of S 

5dP=|\/S = Pj = | VJ G JT, and detP ^ o|, (4.4) 

where P(x^) is any smooth, invertible real matrix function on U. The el¬ 
ements of P have hxed eigenvalues over the coordinate chart W, and P is 
composed of the matrix functions whose ranges are grouped in the similarity 
equivalence classes of the elements of J. Once more, when constructing P 
one must choose the elements of J which exist within the domain defined by 
the determined /dj-parameters, thus, one must respect the validity domains. 
We should remind the reader of the fact that the eigenvalues, ei, and 62 
are the same for the elements of a similarity equivalence class. Referring to 
Eq. fi2.4p we now have 

/ = g^- (4.5) 

In constructing the sets S, and P we treated \/E as a general solution of 
the matrix equation Eq. fl2.16p without considering its relevance to g, and 
/. However, Eq. fl4.5p now brings a constraint on it. Since g, and / are 
symmetric \/S is restricted to the condition 

g^ = E^g. (4.6) 

In other words, not all the elements of S, and P will result in a symmetric 
/ in Eq. fi4.5p when a foreground metric g is specihed. One, has to consider 
the subsets of them whose elements fall into the range of the product of two 
symmetric matrices when squared. For the elements of P, Eq. fl4.6p becomes 

gP-^J^P = {p-^J^Pfg. (4.7) 
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When g is specified, one can choose ten function components of P{x^) freely, 
and determine the other six in terms of these, and the components of g by 
solving the six algebraic equations arising from Eq. fl4.7p . Then, one can 
explicitly construct the symmetric solution / via Eq. fl4.5p in terms of g. 
Thus, we observe that we have P = P[g, x^] indeed. We should remark that, 
although this is the general solution construction method, one may prefer to 
choose simplified, or well-designed forms of P to generate solutions directly. 
For example, if one chooses P = I4, and take a diagonal g one immediately 
obtains diagonal / solutions for the choice of the diagonal elements of We 
now define 

= {{gj) \ f = gX^ \xe S, ^ndgX^ = {X^g}, (4.8) 

and 

= {(gj) I / = giPjf \Pj e P, e^ndgiPjf = {{Pjffg}, (4.9) 

in which, one must construct P[g,x^] in Pj = P~^JP such that it satisfies 
Eq fl4.7p as we discussed above. Next, we will compute the effective cosmo¬ 
logical constant in Eq. fl2.15p which enters into the ^f-equation Eq. fl2.6p as an 
effective presence of the background metric. Taking the trace of Eq. fl2.16p 
leads us to 

/^iCi -|- 2/^262 + 3/^363 = 0, (4-10) 

where we have made use of Eqs. fl2.3p . By using this relation in Eq. fl2.14p . 
and by also referring to Eq. fl2.2p we conclude that 

A = —-(diCi — -(^2^2, (4-11) 

where for a particular solution of the form 

/ = g{p-^[g,x>^]JP[g,x>^]Y = gP-^[g,x^]J^P[g,x>^], (4.12) 

with J & J one must read the appropriate ei(/9j), and 62 (A) for the particular 

choice of J from its parametric structure derived in the previous section. 
This is due to the fact that ci, and 62 are the symmetric polynomials of 
\/S = P~^[g, x^^]JP[g, x^], and under similarity transformations they remain 
the same so they are also the symmetric polynomials of the particular J 
which is chosen to generate the solution. Since, A must stay constant over a 
coordinate chart Li to be able to satisfy Eq. fl2.13p . We will define 

r = r^urp, (4.i3) 
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as the general solution space of Eq. fl2.12p where the set formally corre¬ 
sponds to the elements in the set F^ — F-p with invariant ei, 62 , or invariant 
A over U. We conclude that, the elements of F constitute the decoupling 
solution space of bigravity, with F-p being an analytically expressible subset 
of it that is composed of similarity equivalence classes of matrix functions. 

5 Proportional Backgrounds 

Before we conclude, in this section, we present a class of exact solutions in 
which the two metrics are proportional to each other so that 

ftiu = (5.1) 

These solutions are generated by the elements of Eq. fl4.ip which are pro¬ 
portional to the unit matrix, namely, by {Ui,U 2 ,Vo, K, Xi}. Among such 
solutions which are known as the proportional backgrounds in the literature 
[ 20 ] , as a special class there exists the subset in which both of the metrics are 
diagonal in the same coordinate system. This is due to the fact that in the 
action Eq. fl2.ip both metrics are dehned on the same coordinate patch, and 
in our formalism of deriving these solutions we have kept the same coordi¬ 
nate chart for both of the metrics all through our analysis. After describing 
the general structure of the solutions in Eq. flS.ip we will explicitly con¬ 
struct examples of bidiagonal-metric solutions in which both of the metrics 
are of Friedmann-Lemaitre-Robertson-Walker (FLRW) type. The values of 
the proportionality constants C in Eq. fIS.ip which only depend on the {A}- 
coupling constants for {f/i, f/ 2 , Ro)-^ 1 } are given in the Eqs. fl3.32p . fl3.43p . 
fl3.9ip . respectively. Also, we have C = A = Ci/n, with n = 4, —4 for Ki,K 2 . 
Here, when we computed ei via Eq. fl3.83p for {n,m) = (4,4),(—4,4) for 
Ki, K 2 , respectively we previously found that Ki = K 2 , and we have already 
dehned the matrix K = Ki = K 2 . Thus, one can take C = A = ei[iFi]/4 for 
K. The proportionality constant C corresponding to K also depends only on 
the coefficients {/9j}. As we mentioned in Section two, since the conhgura- 
tion in Eq. fIS.ip satishes Eq. fl2.12p the second of the Bianchi identities in 
Eq. fl2.1ip is directly fulhlled, and the hrst one leads to a contribution of a 
cosmological constant term to the gf-metric equations Eq. fl2.6p . The value 
of this effective cosmological constant can be computed from Eq. fl4.1ip as 

16 

A = -yAC - 4AC^ (5.2) 
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On the other hand, for the solutions of type Eq. fl5.ll) since r = 0 we see via 
Eq. fl2.9p that the /-metric equation Eq. fl2.7p gets no contribution from the 
g-f interaction terms in the action Eq. fl2.ip . Therefore, for the proportional 
background solutions the g, and /-metric equations become 

+ rn'^A)gf,^ = 8nGT^ (5.3) 

and 

f^[Giu - ^A^/^^] = eSvrGTj^^^, (5.4) 

respectively, where and are the corresponding Einstein tensors. 

Since, for the solutions Eq. fl5.ip we have (G-f)^j, = {G^Yj^/C"^ if one chooses 
the fine-tuning of the g, and the /-sector sources as 

A'= ^(A” + m"A), [tI,)\=-YHTI,)\, (5.5) 

then the equations fl5.3p . fl5.4p become equivalent upon the choice of the 
solutions in Eq. fl5.ip . If one of these equations is satished the other is 
also automatically satisfied. The hne-tuning that is introduced above is a 
typical characteristic of the proportional backgrounds I2D1- However, we 
should state that the proportional backgrounds derived here are new and 
different from the already-known ones in the literature. The reader may 
refer to [20] to observe that the known solutions of type fl5.ip are derived 
by excluding the possibility of the vanishing of r which results in a different 
hne-tuning condition that determines the permissible values of C. Therefore, 
we conclude that for any solution of the Einstein equations Eq. fl5.3p which 
are modihed by an effective cosmological constant that is proportional to 
m? with the metric g, matter cosmological constant A®, and the matter 
energy-momentum tensor the metric / chosen as in Eq fl5.ip via the 

introduction of the /-cosmological constant and matter sources as in Eq. 
fl5.5p is a solution of Eq. fl5.4p . Hence, consequently the metrics g, and 
/ = C^g become the solutions of the bigravity action Eq fl2.ip . Next, we will 
construct an explicit example. Let us take the two proportional metrics as 
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FLRW type 


g = —dt^ + 


a?{t) 

1 — 


+ a?{t)r'^d9‘^ + a?{t)r‘^sin^9dip^, 


f 


-C^dt^ + 


1 — 


(ir^ + C‘^a^{t)r‘^d9‘^ + a? {t)r‘^ sin^ dd^p^ ^ 


(|5.6|) 


with any of the proportionality constants C corresponding to the derived 
solntions {t/i, t/ 2 , Vo, K, Xi}. We also take the matter energy-momentnm 
tensor in perfect flnid form 


rpg 

M 


(- 


\ 


-pit) 
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0 

0 


0 

pit) 

0 

0 


0 0 \ 

0 0 

pit) 0 

0 p{t)J 


(5.7) 


where we dehne the matrix Snbstitnting g from Eq. 

and also Eq. 05.71) into Eq. 05.3p leads to the Friedmann eqnations 


and 


/a^2 k SttG a® + m^A 

U = —''-6“^ 


2 a .0x2 k ^ K^ + rn?K 

- = -(-) - — -SnGp ---. 

a 2 


We also have the flnid eqnation 


P = - — {p + p), 
a ' ^ 


(5.8) 

(5.9) 


(5.10) 


for the perfect flnid matter energy-momentnm tensor in Eq. 05.7p coming 
from its conservation law = 0. As nsnal, if one specifies the eqnation 

of state p = wp for the perfect flnid in Eq. 05.7p one can solve the scale factor 
a(t) from the eqnations 05.Sp . 05.9p . 05.1Up . and then determine the metric g 
in Eq. 05.6p explicitly. If fnrthermore, via the fine-tnning that is introdnced 
in Eq. 05.51) we choose 


Tf 

M 


K 

7^ 


(-pit) 
0 
0 

V 0 


0 0 0 \ 

p{t) 0 0 

0 p{t) 0 

0 0 p{t)j 


(5.11) 
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and A'^ = (A® + m?K)/C‘^ then the / metric in Eq. fl5.6p antomatically satis- 
hes Eq. fl5.4p since now this eqnation becomes eqnivalent to Eq. fl5.3p as we 
discnssed above. In this manner, the metrics in Eq. fl5.6p form explicit, and 
exact backgronnd solntions of the bigravity action Eq. fl2.ip for any valne 
of C coefficients of the cases {f/i, f/ 2 , Vq, K, Xi}. We shonld emphasize on an 
important point once more here, before we conclnde. Both in onr example 
of cosmological solntions, and in the more general proportional backgronnd 
cases given in Eq. fIS.ip . the g, and / metrics are dehned in the same coordi¬ 
nate chart. These solntions are similar to the hrst class of solntions classihed 
in [20] . For this reason, we do not have to solve a system of partial differential 
eqnations like the cases in the second class of solntions listed in [20] to hnd 
the Stiickelberg helds which relate the g metric components to the / metric 
ones which are diagonal in different coordinate systems. 

6 Concluding Remarks 

Following the identihcation of a cnbic matrix eqnation which serves as the 
dynamical deconpling solntion ansatz for the two metrics of bigravity, we de¬ 
rived the general solntions of this non-constant-coefficient matrix eqnation. 
We hrst obtained the entire set of solntions of an ordinary cnbic matrix eqna¬ 
tion by classifying them into similarity eqnivalence classes of diagonalizable 
and nondiagonalizable solntions. Later, starting from these cases we have de¬ 
rived the complete set of roots of the actual ansatz equation with coefficients 
also being functions of the elementary symmetric polynomials of the roots. 
Consequently, we obtained the Jordan canonical forms of all the equivalence 
classes of solutions written in terms of the /J-coupling constants of the in¬ 
teraction term of the bigravity action. The decoupling ansatz contributes 
an effective cosmological constant to the foreground metric equation, while 
it leaves the background metric theory as a detached GR. By being able to 
express the elementary symmetric polynomials of the solutions of the cu¬ 
bic matrix ansatz equation in terms of the coupling constants of the theory 
we have also explicitly obtained the effective cosmological constant for each 
class of solutions. Later on, we presented a formal dehnition of the decou¬ 
pling solution moduli of bigravity which forms an important branch of exact 
solutions of the theory. The analytical elements of this effective solution 
space can explicitly be constructed following the solution of six algebraic 
equations. In general, the decoupling solution space tells us which couple of 


31 












the metrics must be chosen to solve the independent Einstein equations in 
two decoupled GR-sectors (one having an always-nonvanishing cosmological 
constant due to the presence of the effective contribution coming from the 
interaction term). Finally, we presented a class of exact solutions in which 
the two metrics are proportional to each other in the same coordinate chart. 
We also explicitly constructed an example of the cosmological proportional 
background solutions of the massive bigravity. We should remark that the 
proportional background solutions obtained in this work are new and they 
differ from the already-existing ones in the literature [20] which are derived 
by excluding the possibility of the vanishing of the matrix r which on the 
contrary sits at the focus of the present work. 

In search for constructing self accelerating cosmologies, the cosmological 
solutions of the ghost-free bigravity have attracted a considerable attention 
in recent years [niC3lElElttHl[l9l[20l[2lll22ll23l^. The known class of 
solutions studied in the corresponding literature can be divided into three 
groups [20]; the class in which both metrics are proportional, a class of 
spherically symmetric solutions in which the background metric is nondi¬ 
agonal, and solutions including diagonal but not proportional metrics. In 
this work, apart from the parametrically-reduced cases, we construct the 
general solution space of the bigravity held equations when the g, and the 
/-sectors completely decouple from each other and simply become GR with 
only an additional contribution of an effective cosmological constant to the 
foreground metric Einstein equations. Beside formally identifying the de¬ 
coupling solution moduli of the theory, we have presented a complete set of 
similarity equivalence classes of analytic solutions up to solving a set of alge¬ 
braic equations. Each of these classes contains functionally inhnite degrees 
of freedom that generate proportional backgrounds, diagonal, and nondiag¬ 
onal solutions. We should also emphasize that, clever choices of ansatz may 
substantially ease the general explicit-solution building method we discussed. 
We believe that, this complete classihcation of the decoupling solutions can 
be a good starting point to study the viable cosmological solutions of bigrav¬ 
ity in a more systematical way. 

A Appendix 

In the Appendix, we present the collection of the matrices constituting the 
set J' dehned in Eq. fl4.ip . Their exact forms whose entries are functions of 
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the {/3j}—coefficients are derived in Section three by substituting the corre¬ 
sponding ansatz into the equation fl3.ip . Firstly, 



Ai 

0 

0 

0\ 


( A 2 

0 

0 

0\ 


/Ai 

0 

0 

0\ 

Ui = 

0 

Ai 

0 

0 

,f/2 = 

0 

A 2 

0 

0 

, 7/3 = 

0 

A 2 

0 

0 

0 

0 

Ai 

0 

0 

0 

A 2 

0 

0 

0 

Ai 

0 


VO 

0 

0 



VO 

0 

0 

A2y 


Vo 

0 

0 

A 2 / 















(A.l) 


where 


with 


Ai = 


Ao — 


P2 + /^sCi — (/^fc^ + 4/32/33C1 + /3| -f 4/3 i/?3)^/^ 

f^2 + /^sCl + (/^fc^ + 4/32/33C1 -|- /3| -t- 4 /?i/33)^/^ 

2^^ ^ 


Cl = 


-6/j2/93T2(9/j|/j|-12/ji/3g)V^ 

3/91 


On the other hand, 


(jD 


(A.3) 


= 


/a; 

0 

0 

0\ 


/A'2 

0 

0 

o\ 


0 

A'l 

0 

0 


0 

a; 

0 

0 


0 

1 

0 

A'l 

0 

1 Uii — 

0 

2 

0 

A'2 

0 

, (A.4) 

Vo 

0 

0 

V 


Vo 

0 

0 

V 



where 


with 


Also, 


A'i = 


A' = 


/^2 + /53C2 — {{ 1^2 + I ^ zC 2 Y + 8 /?i/?3 + 4/32/3302)^'^^ 

W 3 ^ 

f ^2 + /53C2 + ((/32 + hC 2 )‘^ + 8/3i/33 + 4/32/3302)^/^ 

2^^ ^ 


C2 = 


-3/32/33 ±3(/3|/3|-/3i/3|)V2 


/9| 


17 = 


/ —Pi. 

2/33 

0 

0 
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0 

p2 

2/33 

0 

0 


0 

0 

J2_ 

2/33 

0 


0 

0 

0 

3/32 + 2/3361 
2/33 
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(M) 


(A.6) 


(A.7) 
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where ei is a free parameter such that ei G {TZ — {—2(321 ■ 

/O 0 0 0 \ 


U' = 


0 - 


0 

VO 


h 

33 


0 
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0 

^ /33 


0 


0 


0 
0 
0 

2/32+/33e'i 

A 


where e'^ is a free parameter such that e( G {7^ — { 

0 0 0 \ 

0 0 
0 0 


/32 


t/" = 


0 
0 
V 0 


/52 

0 

0 


(A.8) 

—‘i(32/(3z}'\- In addition, 

(A.9) 


For the Id-series we have, 
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where 


R = 




2/^3' 


The X-series elements in are 
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If we turn our attention on the X-series we have 
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